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Optimal regularity and nondegeneracy of a free 
boundary problem related to the fractional 

Laplacian 

Ray Yang 
Abstract 



We discuss the optimal regularity and nondegeneracy of a free boundary problem 
P_i ' related to the fractional Laplacian. This work is related to, but addresses a different 

. problem from, recent work of Caffarelli, Roquejoffre, and Sire [5J. A variant of the 

(-H I boundary Harnack inequality is also proved, where it is no longer required that the 

■ function be along the boundary. 

1 Introduction 

In this paper, we prove certain local properties (regularity and nondegeneracy) near the free 
boundary for minimizers of an energy functional related to the fractional Laplace operator 
CS| ' with a penalty term. This work is inspired by recent work Q on a similar problem, and we 

explain the broader background first. 

In the theory of one-phase free boundary problems arising from the minimization of 
energy for the classical Laplacian, Alt and Caffarelli [2J analyzed minimizers of the energy 
7(a) = J (Va)^ + Xu>odx subject to non-negative Dirichlet data, while the study of the free 

^ ! boundary arising from minimizers of the energy J(u) = J ^^j- + udx with non-negative 

^ \ Dirichlet data is encompassed by the study of the obstacle problem. An intermediate case 

is the case studied by Alt and Phillips [3, which is that of the free boundary for minimizers 
of J{u) = J + u^dx, where < y < 1 . In a heuristic sense, we can view the Alt- 
Caffarelli problem as the case of y = 0, and the case y = 1 as a special case of the obstacle 
problem. 

The problem we study is the analogue of the problem of Alt and Phillips for the frac- 
tional, rather than standard, Laplace operator. The current article only covers the regularity 
and nondegeneracy of energy minimizers, and is thus properly the analogue of Phillips' 
work in [17J and part of [16J. 

The fractional Laplace is a nonlocal integral operator, taking the form 



(-A) u(x) = -^^—dz. 

Jr" \x — z\ 
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This operator has a corresponding energy given by 

r (u(x)-u(z))\ , 

J{U) = I — ; — dZClX 

JR"xR" z\ 

This latter term is a nonlocal energy, and not very easy to manipulate. In (7], CafFarelli and 
Silvestre introduced the notion of extension to one extra spatial dimension and examining 
a particular PDE on the upper half-space, with the fractional Laplacian being equivalent to 
the Dirichlet-to-Neumann map at the boundary. To be precise, consider (x, y) e R" x 
and a function u on the upper half-space satisfying 

v(y-^''vu) = 

Then 

(-ATu{x,0) = -Cn,Aimy^'^'^dyU(x,y) 

The natural energy corresponding to the second-order equation on the half-space is then 
7(m) = J y^~^"'\V u{x, y)\^dxdy . This is an energy where it is easier to study the purely local 
properties of its minimizers. 

The extension characterization of the fractional Laplacian has been used to study both 
the obstacle problem ([9J) and the equivalent y = case (L5J). We will apply it to study the 
intermediary case, that is to say, minimizers of energies 

J(u)= f y^-^'^\Vu(x,y)\'^dxdy + [ u^dx 

J J{y=0] 

in subsets of the upper half-space with parts of their boundary lying along 3^ = 0, where we 
have < y < 1 . Since this is a study of the one-phase problem, we assume non-negative 
Dirichlet boundary conditions. The second term in the energy penalizes non-zero values of 
u along the hyperplane {y = 0}. Hence, we can consider separately the zero set of u (called 
the contact set), and its positivity set. Restricted to {y = 0}, the interface between the two 
is the free boundary. 

We prove that in any neighborhood of a free boundary point, our energy minimizer u 
lies in the class with a Holder seminorm that depends only on the distance to the free 
boundary, where /3 = ^ is the critical scaling exponent for the problem. This is called 
the optimal regularity of u, since we also prove the non-degenerate nature of u, namely, 
that in any ball of radius r about a free boundary point, sup^ u > Cr^ for a constant C that 
depends only on n, cr, and y. 

In the course of proving optimal regularity, we prove and use an improvement on the 
boundary Harnack inequality of CafFarelli, Fabes, Mortola, and Salsa BUl, which may be of 
interest even to those not working in free boundaries. 

The ideas behind the proof of the optimal regularity of energy minimizers for fractional- 
order cases can be extended to a proof of optimal regularity for the second-order case, 
which was first proved by Phillips [fTTl . Since the proof for the second-order case illustrates 
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the ideas in a less involved setting than the fractional-order case, we provide it as well. The 
key ingredient for optimal regularity is the construction of a lower barrier, or subsolution, 
for the energy minimizer which is strictly positive at the center of a ball when the values 
near-by are "too large," thus, for a free boundary point to exist, the growth cannot be too 
great. 

2 Preliminary considerations 

In this section we identify some technical points of interest. First, we prove that minimizers 
of the energy exist. Second, we identify the scaling associated with the problem. Third, we 
list certain properties of the equation and minimizer that are known and will prove useful 
to our analysis. 

2.1 Existence considerations and some definitions 

We consider, in B+ = Bi Ci {y > 0}, minimizers of the energy 



negative Dirichlet conditions on dBi r\ {y > 0}, where a = I - 2cr, < cr < l,0<y< 1. 
For sake of convenience, we denote T = Bi f) {y = 0}. 

The energy can be interpreted as an averaging term which "lifts" the solution towards 
the boundary conditions, and a term which punishes u for being nonzero at y = 0, causing 
it to "stick." The set {u = 0}, which necessarily lies in {y = 0}, is called the contact set of u. 
The interface between {u = 0} D {y = 0} and {u > 0} Ci {y = 0} is called the free boundary. 

Existence of minimizers is assured by the usual methods: consider a minimizing se- 
quence for the energy. The first term of the energy is lower semicontinuous with respect to 
the norm for the usual reasons. The second term is continuous with respect to the norm for 
L^(T) From the extension result of Caffarelli-Silvestre we know that the trace of functions 
lying in H^{B+, a) lie in H"'(Y) whence we apply the usual Sobolev embedding of 
inside L^. 

We will use X — (x, y), where x 6 M." and y g R+. 
2.2 Scaling of the problem 

Put briefly, we seek the scaling that preserves minimizers of the energy 




in the space 



(B+, a), with seminorm given by \\u 



jg y"(\Vuf')dxdy. We impose non- 
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We consider the scaling 

1 



wix,y) = —uiAx,Ay) 



We find tliat 



7(w) = l r y"A^-^'^\Vu{Ax,Ay)\^dxdy + f A'^^u^dx 
2 Jn Jr 

which, after the change of variable, scales to 

^-a-n^2-2ii]_ C y'\^uUxdy + A'^^'"^' r u'dx 
2 Jin Jay 

Setting the exponents equal, we find that 

la 



2.3 Some other useful properties 



2-r 



The Euler-Lagrange equations for J{u) tell us that, in a distributional sense, the minimizer 
u satisfies 

V • (y^Vu) = 

in the interior of 5+, and 

Mmy'^dyU = yu^~^ 

along r wherever u > 0. Caffarelli and Silvestre [|7l showed that, in the upper half space 
R"'^^ the eff"ective normal derivative operator is equivalent to the fractional Laplacian of 
order cr: 

Vimy"dyU = -C(-ATu 

Fabes, Jerison, Kenig, and Serapioni ([[T2 l . [fT0l .[|TT | ') extended the De Giorgi-Nash- 
Moser theory of divergence-form elliptic equations to degenerate elliptic equations with 
Muckehnhoupt A2 weights; these are equations of the form 

V • (A(x)Vm) = 

where the matrix A satisfies 

Aw(xm' < J]a,j(x)^i^j < Awixm'' 
for ^ e R", with weight functions w satisfying 



— r w(jc)<ijci( — r — 

\B\Jb I\\B\J w(. 



dx\ < C 

for all balls B. In particular, such properties as the strong maximum principle. Holder 
regularity, and the Harnack inequality all hold. There are certain other properties, such as 
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the De Giorgi Oscillation Lemma and a more specific form of the De Giorgi-Nash-Moser 
Hamack inequality, which follow directly from their work but are not explicitly stated. A 
discussion of those properties has been provided in iB]for the reader's convenience. 

Finally, the generalized Hopf lemma for cr-harmonic functions in R" (first stated in f5^) 
is useful: 

Lemma 2.1. If a smooth function v{x) satisfies (-A)°"v = in some smooth domain Q c R", 
V being non-negative and positive in the interior ofQ, and if there exists a point xq e da 
where v(x) = 0, then there exists C such that v{x) > C{{x - Xq) • y(xo))°" where v{xq) is the 
interior normal to dQ. at xq. 



3 Optimal regularity for the 2nd order case 

The optimal regularity for the problem in the 2nd order case was first obtained by Phillips 
[[TTll . Our method for proving the optimal regularity of the fractional case can be adapted 
to give an alternative proof for the 2nd order case. The main intuition behind our proof in 
the 2nd order case is free of certain technical issues that occur in the fractional case, and so 
we present it here first. 

We obtain the optimal regularity of the energy minimizer u to the energy functional 

r ivmP 

J(u) = I — ^ — I- wdx, 

showing that u 6 C^'^"^ where p = the scaling factor obtained by a calculation like that 
in ^2.2[ As with the fractional case, we assume the boundary data is non-negative, which 
allows us to assume the same for u. Notice that the Euler-Lagrange equations tells us that, 
when u > 0, u satisfies 

Au = yu^~^ 

As with Phillips, we seek to prove 

Lemma 3.1. There exists a constant c^in) such that if 

udS > cor^ 

then u > inside B^. 

Since the scaling u^ix) = jpu(Ax) preserves minimizers, we need only show this for 
r = 1, and scaling would take care of the rest. 

Our proof works by showing that, when the average on the boundary is sufficiently 
large, a subsolution, or lower barrier to the energy minimizer, can be constructed which is 
wholly positive in the interior of B^. There are two main stages to the proof: first, we detail 
what it means to be a subsolution, second, we construct a subsolution with the desired 
properties. 
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3.1 Subsolutions 

We say that a function w is a subsolution, or lower barrier, to the energy minimizer, if w 
satisfies 

Aw > Mw^~^ 

whenever w > inside Br, setting w = u along dB^ for Dirichlet boundary conditions, with 
M a large constant to be determined later. This terminology is natural because, as we shall 

see, u>w inside B^. 

Let V = max(M, w), and consider the difference of energies given by 

J(u) - J(v) = J\ ((^")^ - (Vv)^) + - v^dx 
Since u is the energy minimizer, we require that J(u) - J(v) < 0. However, we know that 

i J (Vuf - {Vvfdx = ^Jv(u + v)- V(m - v)dx 

1 r 

^ ^ f (v - u)(ju^-^ + Mv^-^)dx + ^ f Mv^dx 

^ Jv>u>0 ^ Jv>u=Q 



We compare 
with 



(P{s) = \{s-t){yr'+Ms-y-') 



for s > ? > 0. Clearly, 0(0 = il/{t) = 0, and a bit of calculation assures us that (pXs) > (//'(s) 
for all s if we set M >2. Thus, 

m-m = /i((v„,^-(Vv)V„'-vr.. 

^ ^ r (v - u)(yu^-^ + Mv'^-^)dx + f Mv^dx 

^ Jv>u>0 ^ Jv>u=Q 

> 

with equality holding in the last statement only if v = m. Hence, u>w, with M = 2. 

3.2 Construction of a positive subsolution 

Our goal is to create a positive function w on the unit ball, such that Aw > 2w^~^ We will 
define w in three parts: 

w(x) = wi(jc) + W2(x) + w^(x) 
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Let ri(x) be a radial, non-negative C°° function satisfying 77 < 1 everywhere, 77 = 1 when 
> i> ^7 = when |x| < ^, and IA77I < C and IV77I < C for some constant C. We define 

w,(x) = A(r](x)(l-\x\f + (l-r](x))) 

and we claim that when |x| > |, we have 

for the correct choice of A. 

In the region in question, it is easy to see that 

AwiW =/3(J3- 1)(1 - r/-2 - (n - 1)^6 ~ ""^ ' 

r 

We take the ratio of Awi with wi in the region under concern, and we see that 



whence it is clear that a sufficiently large value of A will suffice. 
We set 

wiix) =MW^ - 1) 

where we pick jji suflaciently large so that Awi > -Awi + 1 everywhere inside Bi_. It is clear 
from our design that -Awi is bounded inside the region in question. 

Finally, we let W3 be the function which is harmonic inside Bi, with the same boundary 
values as the minimizer u along dBi. 

We claim that when j^^^ udS is sufficiently large 

1. 1 > 2w^-i on 57. 

8 

2. w > everywhere on Bi 

The Hamack inequality tells us that, on 5?, we have 

wj,{x) >C-r udS 

JdBi 



To prove the first, it suffices if -^^ udS is so large compared to // so that 2/in > (C udS y~^. 
To prove the second, we bound W3 from below by a suitably scaled truncated fundamental 
solution, and then 

(f)"-2- 1 \\X\" 



when f udS is sufficiently large. 



Hence, we have w > everywhere, and on 5?, we have 

Aw > 1 > lw^~^ 

while on 5 1 \Bi we have 

Aw > Awi > 2w\'^ > 2w^~^ 

and we are done. 

4 Optimal regularity for the fractional case 

The goal of this section is to obtain the optimal regularity of energy minimizers u. In 
particular, we seek to show that u grows away from the free boundary like a power of the 
distance. To be precise, u{X) < Cd^, where d is the distance of X to the free boundary, and 
jS = ^ is the scaling factor obtained in ^2.2[ 

As corollaries, we obtain some regularity results: restricted to {y = 0}, u lies in the 
Holder space C^. If > 1, we will prove that u e C^'^"\ which by abuse of notation we 
will still refer to as C^. In the interior domain where j > 0, we still obtain u e when 
< 1 . When fi> 1 , we find that u e C" for any a < \. 

To obtain optimal growth, we consider a point po which is at some distance (normalized 
to 1) from the nearest free boundary point, which we will take to be 0. We will use a variant 
of the boundary Hamack inequality due to Caffarelli to compare the value of u{po) with 
some point pi in the interior of 5+, specifically, showing that u(pi) > Mu(po). We can 
then use the regular Harnack inequality in the interior to show that, in a smaller ball about 
the free boundary point, the boundary values are controlled by u(pi), and hence by u(po). 
We then prove that, if the boundary values in the upper half ball are too large at the right 
scale, then m(0) is strictly positive, meaning that u(po) cannot be too large. Subsequently, 
rescaling obtains the desired regularity. 

Our main tools to prove optimal growth are a variant of the Boundary Hamack In- 
equalit}0, and a lemma stating that if a particular weighted integral along the boundary of 
the a half-sphere is sufficiently large, then the minimizer of the energy taking boundary 
conditions along the sphere has a positive value at the center. 

Theorem 4.1 (Variant Boundary Hamack Inequality). Let u be a non-negative solution of 
the equation V • (y"Vu) = in B^, which satisfies \imy^Qy"dyU > along {y = 0}, taking on 
some continuous boundary values along {y = 0], with u{0, ^) = 1. Then inside we have 

2 

u < M for some constant M{n, cr). 

^There are two types of results which are, confusingly, both called the Boundary Harnack InequaUty in 
the literature. In addition to the result here, which states that values in the neighborhood of the boundary are 
uniformly bounded in terms of the value at an interior point, there is a closely related result which states that 
for two solutions which are both along a stretch of the boundary, their ratios are locally Holder-continuous. 
We follow the naming convention of Caffarelli and Salsa |8| and call the first result the Boundary Harnack 
Inequahty, and the second result the Boundary Comparison Principle. 
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Lemma 4.2 (Minimizers with large averages are positive at the center). Let u be a mini- 
mizer of the energy J{u) inside Br C\{y > 0}, taking non-negative boundary values along 
dBr n > 0}. 3co > such that, if 

u\dBr\\y>^_] ^ Co 

then we have 

u(x,y) > Oy(x,y) e 5^ n {y > 0} 
and in fact there exists a constant c such that 

u(x,y)>c^(x,y)eBrr\{y>0} 

Together, these suffice to prove our resuk, namely: 

Theorem 4.3. There exists a constant K such that in any B^Xo) n F, where Xo is a point 
such that u{xq) = 0, such that 

\u{x) - u{xq)\ < K\x - Xof 

where /3 = 

Proof. Without loss of generality, let be a point such that a(0) = 0, and X* be a point 
such that 1X1 = 1. We claim that < K. Suppose this is not true, that is to say, we can 
make as large as we wish. Then by the variant boundary Harnack inequality applied 
to 52(0) we have 

u(0, 1) > 

where M is the constant from the variant boundary Harnack inequality. By applying the 
DeGiorgi-Nash-Moser Harnack inequality to u about we discover that we have 

u{x,y) > Cu(Xt)y 



whence, by invoking Lemma 1421 we have m(0) > 0, a contradiction on our original as- 
sumption. Thus, there exists a constant K such that m(X,) < K, as desired. By rescaling the 
problem, we recover our desired result. □ 

Corollary 4.4. Let u be an energy minimizer in a subset o/R""^' containing B"^, with a 
free boundary point. Then considered as a function along the set [y = 0}, u is a function, 
with \\u\\cfi(Bi) ^ C, where C depends only on cr, y, and n. 
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Corollary 4.5. Let u be an energy minimizer in a subset o/R.""^ containing 5| with 
a free boundary point. Then in B\, u is a function, with \\u\\cP(b+) < C, where C depends 

only on cr, y, and n, if p < 1. Ifp <\,u is a function for any a < 1, with 



|m|Ic-(b+ ^ C{a,y,n,a) 

1 



The proof of these statements, and a discussion of the norm estimates, are covered 
in 1431 
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4.1 Variant Boundary Hamack Inequality 



The proof of this variant of the boundary Hamack inequality follows the same lines as 
the standard proof of the boundary Hamack inequality provided by Caffarelli et alia 
The proof uses two classical facts from the De Giorgi-Nash-Moser theory, which was ex- 
tended to the theory of degenerate elliptic equations with A2 weights by Fabes, Kenig, and 
Serapioni [12], a class that includes the equation V • (y"Vv) = 0. 

The first fact is the De Giorgi-Nash-Moser Harnack inequality, which states that for a 
non-negative solution inBi, 

sup u < c(\ - ry inf u 

Br 

where p > 0. 

The second fact is the De Giorgi oscillation lemma, which says that a subsolution v in 
Bi satisfying 

• V < 1 

• |{v < 0}| = a > 

has the property that 

sup V < fi{a) < 1 

Bl/2 

We proceed by contradiction. Let m(0, ^) = 1, and extend u over the line j = by even 
reflection. Suppose there is no M which can bound values of u inside the half -ball. Then u 
achieves its maximum in it Mq > M, at some point Xq = ixo,yo). The Hamack inequality 
tells us that the distance to the boundary, yo, satisfies 

3^0 < ^0 = 

We now proceed with a construction we repeat for each successive value of n, starting with 
n = 0: 

Consider now BKd,X^n, 0) (the projection of X„ to the plane y = 0), for K large, greater 
than, say, 4. For points satisfying y > 2dn, we have, by the Hamack inequality, that 

M 

u(X) < c{2d„r = 

The set {y > 2d„} has measure at least a fixed fraction of BKd,X^n,0), independent of K. 
Thus, if we let M„+i = sup^^^ (-^.^ q) we see that inside, say, B2dX-Xn,0), by the oscillation 
lemma, we have that 

M„{l-2'P)< sup (u-^)<ij(K)(Mn^i-^) 
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with fi —> as K becomes large. Thus 

M„+i > Mn(2-P + ^—^) 

We pick K sufficiently large that the factor on the right hand side is some fixed positive 
A> I. We let Xn+i be the point where = M„+i inside BjcdX^n, 0). 

Thus we have a sequence of points X„. Notice that K does not change, and hence neither 
does A. h%n ^ oo, we have 

u{Xn) > A"Mo oo 

while 

y„<dn = (cMoyh'r 0. 
The distances between the points satisfy 

\Xn+l - Xn\ < Kdn 



and so the sequence has 

d{X^, X„) <Ky d„ < if ^£^1^ 
which can be made to converge inside B9_ for M sufficiently large, giving us a sequence of 

16 

points X„, with limit points where u blows up along y = 0. This contradicts our original 
assumption that u continuously assumes values along the boundary {^^ = 0}. 



4.2 The center is positive when the boundary is large 

Lemma l4r2l consists of demonstrating that when co is sufficiently large, a subsolution which 
is purely positive in can be built, which serves as a lower barrier to the solution. 

4.2.1 Conditions for a subsolution 

We seek sufficient conditions for a function to be a subsolution of our variational problem. 
One way to do this is to show that, for a subsolution w, where m > w, we can improve the 
energy: if v = max(M, w), then 

J{u) - J(v) = ^J / (|VmP - |Vvp) dydx + - ^^^^ ^ 

since u is the energy minimizer. Clearly, the second term is negative; our approach lies on 
setting conditions so that the first term dominates the second. 



11 



We assume u, v sharing the same Dirichlet boundary conditions along dB, and integrate 
by parts: 

J ^/(IVwp - \Vv\^)dxdy = \J /(^w + Vv) • (Vm - Vv)dxdy 

= ~l: A" ~ ■ (/^(" + v))dxdy - \ [{u - v) limy^dyiu + v)dx 
2 J 2 Jr .v^o 

We define v = max(M, w), where w satisfies 



\imy"dyW > Mw'^'^ 

y-,0 ^ 



along r, and V • iy^^Vw) = in 5i . Then lim^.^jo dyV ^ Mv^ ' on those portions where v > u, 
whence we can write 

\ y\\Vu\^ -\Vv\^)dxdy> \ (v - u){yu^-^ + Mv^-^)dx + f Mv^dx 

J Jrn{v>M>oi Jrn{v>u=o) 

Recall now that 

7(m) - J(v) = ^ r /(|VmP - \Vv\^)dxdy + [ - v^dx 
2 Jb Jr 

Since u is the energy minimizer, we need for this term to be negative. 
We consider the functions 

iff(s) = s^-f 

and 

Clearly, 0(0 = tff{t) = 0. We now examine their behavior in the range < t < s. When 
s > t, 

i]j'{s) = ys^-^ < (p'(s) = ^yf-^ + yMs^-^ + y(l - y)Mts^-^ 
. Thus, (f){s) > ifr{s) when s > t, and we can write 

J{u)-J{v) = ]: \ yWVuf -\'^v\^)dxdy+ \ u^-v^dx 
2 Jb Jr 

> -[ (v - u)iyu'^-^ + Mv^-^)dx + I- f Mv'^dx+ fu^-v^dx 
2 Jrn{v>«>0) 2 Jrn|v>n=0| Jr 

> 

with the last equality being strict if v differs from w on a set with positive measure. This is 
satisfied if we set M = 2. 

Hence, whenever such a w exists, we can decrease the energy of u, a contradiction on 
the definition of u as the energy minimizer. The construction of such a w, which is positive 
on 5i and greater than a fixed constant on 5i , is thus sufficient to prove Lemma |42l 

3 6 
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4.2.2 Construction of such a subsolution 

We want our subsolution w to have three properties: we would like our w to take the same 
values as u along dB Ci {y > 0}, we would like it to satisfy the conditions 

V • (/Vw) = 

in B+, and along Y we would like 

limv"(9,.w > 2w^~^ 

wherever w > 0, and finally we want w > O'mB^ n F. We will define our w in two parts. 

W = Wi +W2 

We set wi by setting, for x e R", 



Let 



(0 \x\ > i 

l-(l -3|jc|/-2^ \x\<\ 



-c f 



be the Riesz potential of ifr. We state a technical lemma relying on classical results in the 
theory of fractional integration and Riesz potentials, which leave to appendix lAl 

Lemma 4.6. (/2o-iA)(-'c) is well defined and continuous as a function, radial, has fractional 
Laplacian equal to i//(x), and furthermore, there exists 6 > such that 

|(/2.(A)(r) - (/2.iA)(^)l < C(l - 3rr 

for ^>r = \x\>^-6 where min(yS, 1) > a > cr. 

We let b(x) be equal to /20-iA on R" \Bi, and have (-A^b = inside Bi, and then we 

set 

w(x) = (ho-t/^Xx) - b{x) 

b is the solution to the standard Dirichlet problem for the fractional Laplacian; its ex- 
istence is guaranteed by the standard theory (see, e.g., Landkof iflSlD . Notice that w is 
cr- subharmonic. This means it is negative inside B^, and outside of it. Furthermore, 
the maximum principle for cr-harmonic functions (Lemma [2TI) applied to b{x) tells us that 
there is a constant such that 

\b{x)-b{^)\<C{\-3rr 
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Now we let 

,,2o- 



wi(x,j) = C„,^ —dz 

J iix-zy+y^)^ 



where z ranges over R.". This is, of course, the Poisson kernel for the fractional Laplacian 
convolved with w, giving us a wi that satisfies V • (j"Vwi) = in the interior, which takes 
on the values of w along {y = 0}, satisfying \imy^Qy"dyWi{x,0) = if/(x) (by the extension 
result of CafFarelli and Silvestre fl\). 

For the sake of future estimations, it is helpful to bound -wi from above by an auxiliary 
function. We let qo = 2 sup(-wi), and we let 

\qo \x\ < ^-6 



\x\ > \ 



q(x) = 

and let q satisfy (-A)""^ = on the annular ring ^ - 6 < \x\ < ^. The comparison principle 
for fractional-harmonic functions then tells us that q > -wi on Y. We extend q to R."^^ in 
the usual way via the Poisson kernel. 

Q(x,y) = Cn,<T —dz 

J ((x-zy+y^)^ 

Thus, we have 

Proposition 4.7. Q(x,y) > \wi(x,y)\ in the upper half-space W^^. 
We set W2 with boundary conditions 



\u-WiiX) XedBif\[y>Q} 
lO XeT\Bi 



W2{X) = 

and let it satisfy the problem 

V • (y''Vw2) = 

when X e B+, and 

\imy"dyW2 = 

whenX 6 rn5i. 

3 

Now we need to estimate properties of w = wi + W2, which we do by comparing W2 to 

Q 

Proposition 4.8. For any A > 0, a sufficiently large value of cq will ensure that W2{X) > 
AQ{X) inside 5+ 
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Proof. Since both functions satisfy V ■ y"(yv) = inside 5+, it suffices to examine their 
relative behavior along C\{y > 0} and along Y. 

Along n > 0}, the boundary comparison principle (see IfTTTl for a proof in the case 
of A2 weighted degenerate elliptic equations) tells us that 

u{x,y) > Cco/'^ 

since a > Cq when y > \- We also have from the formula that that Q{x, y) < Cq^y^"' . Hence, 
we just need co to be sufficiently large. 

The behavior along F is a touch trickier. We divide our analysis of the behavior of W2 
along r into two parts: the first part concerns the interior of B^_g, where 6 is from Lemma 
14. 6[ and the other in the thin annular ring | - 5 < |x| < |. Clearly, F \ 5^ is not taken care 
of, since W2 and Q are identically there. 

We note that W2\dBr\{y>o] > 0, so that we can apply the Harnack inequality in the interior. 
We bound W2 from below by a function W2, which we define as follows: let 



W2{X) 

and let it satisfy the problem 



\u-Wi{X) X e^B^(^{y>Q} 
[0 XeY 

V ■ (/Vvv2) = 



Clearly < VV2 < W2 in the domain. Since we know W2 in {y > \} is greater than cq, it 
follows that so too W2 at interior points, such as, say, X = (0, |), is linear in cq, and hence 
so is W2- We apply the Harnack inequality to W2 inside the ball 5i to see that W2 > Ccq in 
Bi_g can be made as large as we wish, where 6 is from lemma 1461 Thus, inside Bi g Pi F, 
we can choose cq so that W2 > Q. 

In the annular ring proper, both Q(x,y) and W2(x,y) satisfy limy^o j'^^^yV = 0, whence 
we can invoke the Hopf lemma to see that W2 > Q. □ 

Corollary 4.9. 

w = wi + W2 > (A - 1)2 
and hence by making A sufficiently large we can make w > Cqo in Bi. 

6 

We close our construction with a lemma, which shows that w has all the desired prop- 
erties of a subsolution. 

Lemma 4.10. 

wi + W2 > 

and for A from the previous proposition sufficiently large ( which is really to say for cq 
sufficiently large), we have 

limy^^vW > 2w^~^ 

along Y, wherever w 0. 
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Proof. Since Wi>-Q and W2> AQ, we have 

W = Wi + W2 > W2 - 2 > 



On r n 5i c, we have 

3 " 



Y\my"dyW = Y\my"dyWi 



By setting X sufficiently large, we can attain 



On the annular ring, we invoke Lemma [27T] to see that there is a constant c such that 

q{x) > c(l - 3|jc|r 

whence we derive the relation 

w{x) >{A- l)c(l - 3\x\f 
By setting A sufficiently large, we can make 

2FT(i - i)c(l - 3|jc|r > (1 - 3|x|/ 
in the entire annular ring, whence we can attain 

limy^^vW = limy^^vWi 

= (l-3W)^-2- 
= (l-3Wf^^') 

> 2((^-i)c(i-3wrr' 



□ 
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4.3 estimates for u 



The goal of this subsection is to provide a proof for Corollaries 14.41 and 14.51 

For the first, we will do this by analyzing the effective equation satisfied by u, restricted 

to r, in the neighborhood of a free boundary point. The estimates follow the spirit of the 

analysis conducted in Section III of [17J: we will first show that appropriate Holder norms 

of u satisfy certain pointwise estimates in terms of the value of u itself, and then put these 

estimates together to obtain a uniform estimate. 

For the second, we will follow a similar procedure, first using interior estimates to get 

pointwise bounds on Vm when y > 0, and then tie these together with the estimate along 

{y = 0) to get a uniform estimate. 



4.3.1 Along {y = 0} 

Lemma 4.11. {m > 0} n F is open with respect to R", on which u satisfies 

Yimy''dyU = yu^~^ 

and furthermore u e C°°{{u > 0}), such that the tangential derivatives of u, which we 
represent by V^u, satisfy 

\VMpM<C(u(p,0))'-^ 

and, moreover, that the tangential second derivatives of u, which we represent by V^^u, 
satisfy 

\V,xu(p,0)\<Ciuip,0))''^ 
where p is any point in [u > 0} n {y = 0]. 

Proof. If p e F is some point such that u{p,0) > 0, then the variant Boundary Harnack 
inequality tells us that u{p, 6) > Cu{p, 0). If we make the usual dilation by A about p, we 
see that u^(x,y) = jpU(A(x - p) + p. Ay) satisfies 

Ua(.p, 1) > c^uip,0) 

For A sufficiently small, cA~^u(p, 0) can be made larger than the constant needed in Lemma 
14. 2[ whence > C > in 0), or, in the original u, we can say u > CA^ > in 

6 

Ba(p,0) for A sufficiently small. 

6 

Hence, the set m > is open with respect to F, and, on every set D compactly contained 
within {u > 0}, is bounded away from zero, hence 

\imy"dyu = yu^'^ 6 V"{D) 

Consider now B = Bi{p,0). We let Ui be the Riesz potential of -JuV^Xb, and set 
U2 = Ux - Ui. Then on B we have 



Yimy'^dyUi = yu\ 



r-i 



17 



while 

\imy"dyU2 = 

and both of the m, satisfy V • (y^Vuj) = 0. 

Since the tangential derivatives of U2 also satisfy the same equations as U2, we have that 
U2 e C~(5), with 

|V,M2(aO)|<C 

and 

\^xxU2(p,0)\<C 

in Bi{p,0), by the estimates found in Similarly, we can use the potential-theoretic 
estimates found in filSil iteratively to show that ui e C°°{B), with 

and 

|V,-,-i^i(AO)|<C. 

Hence, after rescaling, we can say that, for the tangential derivatives of u, we have 

|V,M(p,0)|<C/-^ 

and 

How small need A be? Our condition was that cA~^u{p, 0) > cq, whence we see that 
A = (Cu(p,0))i^ suffices. The conclusion follows. □ 

Up to now, it has been possible to treat the cases where /3 > I and < 1 as if they were 
the same. For the remaining two theorems, we have to recognize the difference. The result 
here is proved very much the style of [fTTl and ||6l. 

Theorem 4.12. Suppose /3 < I. Then there exists a K = K{6,n,j3), such that if X\,X2 e R" 
are in a 6-neighborhood of the free boundary, we have 

\u{XuO) - U{X2,0)\ < K\Xi - X2f 

If li > 1, there exists a K = K(6,n,fi), such that if X\,X2 g R" are in a 6-neighborhood of 
the free boundary, we have 

|V,m(xi,0) - V,M(X2,0)| < K\x, - X2t' 

In either case, since away from the 5- neighborhood of the free boundary, u e C°°, this 
means we can put the two together to get a uniform norm for u. 
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Proof. As in the previous lemma, we notice that there is a constant C\, such that if u{xi , 0) > 
Ci, then the variant boundary Harnack inequality tells us that u satisfies the conditions for 
Lemma 1521 and hence u > C2 inside Bi{xi,0). Rescaling this statement, we have that if 

6 

u(xuO) > Cif^, then u > C2T^ inside Bl{xi,0). 
We now consider three cases: 

1. uixuO) > Ci(6|xi - X2\f and \xi - JC2I < f 

2. u(xi,0) > Ci(6Ui - xjlf and \xi - jcjI > f 

3. u{xi,0),u(x2,0) < Ci{6\xi - X2\'f 

i j Consider the line segment joining xi and ^2, with r = 6|jci-jc2|- Since u{xi,0) > C\r^, 
we have u{x,Q) > €21^ inside B^{xi,0), which happily is precisely B\x^-jc-,\(xi,0). Hence, 
when /3 < 1, the mean value theorem applied along this line segment tells us 

\u(xi,0) - u(x2,0)\ < \V^u(x',0)\\xi - X2\ 

where x' is some point along our line segment. By applying the estimates from Lemma 
14.1 1[ we have 

|V;,m(a:',0)||jci - X2I < C(m(/, 0))^|jci - X2I < Qxi - X2f 
When /3 > 1, we consider instead 

\VMxi,0) - VMx2,0)\ < \V,Mx',0)\\xi - X2\ 
and by applying the estimates on the tangential second derivatives, we have 

\V,Mx',0)\\Xi - X2\ < C{U{X' ,Q))^\Xy - X2\ < C\x^ - ^2^^ 

2 j In this case, we simply say directly that, if < 1 , 

\u{xi,Q) - u{x2,G)\ < \u{xi,0)\ + \u{x2,0)\ 

< < C\Xi - X2f 

where we invoke Theorem 14. 31 on the last step. If jS > 1, we say 

|V^m(xi,0)- V^m(;c2,0)| < |V^m(xi,0)| + |V^m(x2,0)| 

< C(u{xx , 0))^ + C{u{x2, 0))^ 

where we invoke Theorem 14. 3 1 on the last step. 

3) The calculations are exactly like case 2, only instead of invoking Theorem 14. 3 1 to 
bound u pointwise, we invoke the hypothesis. □ 

From this result. Corollary 14.41 is obvious. 
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4.3.2 The estimates when y > 

Note first that inside y > 0, u e C°°, since V ■ {y"Vu) = is uniformly elliptic with smooth 
coefficients on any compact subset contained within {y > 0} (with differing ellipticities, of 
course). We can thus assume that u is smooth far away, and concentrate on its behavior for 
small values of y. 

We start with an elementary lemma that gives us pointwise estimates on the derivatives 
of u via rescaling: 

Lemma 4.13. Let ube a non-negative function satisfying V ■ {y'^Vu) = inside 5r n {y > 0} 
for some large R . Then there is a constant C depending only on n and cr, such that 

C 

\Vu(xo,yo)\ < —u(xo,yo) 

yo 



and 

C 

yl 

whenever yo > and Byo{{xQ,yQ)) c Br. 



|D^M(xo,yo)I < ^u{xo,yo) 



Proof. Suppose first that = 1 and Bi{xQ, 1) is inside Br. Then inside this ball, y'^ is a 
bounded, C°° coefficient, so the standard regularity theory for weak solutions gives us the 
estimates 

|Vm(;co, 1)1 ^ Cu{xq, 1) 



and 



\D^u{x^, 1)1 < Cm(xo, 1). 



For general y, we simply consider the rescaling w{x,y) = u{xq + {x - xo)yo,yQy) and write 
the estimate for w in terms of w. □ 

Next, we provide a boundary estimate on the growth of u away from the line y = 0. We 
choose nice constants for the varius radii and the lines, bearing in mind that we can rescale. 

Lemma 4.14. Let u be an energy minimizer inside B^r\{y > 0} with nontrivial free bound- 
ary. Then there exists a constant C such that, for y < I and (x, y) e Bt„ we have 

\u{x,y)-u(x,0)\ < Cy^ 

Proof. If u{x, 0) = for any (x, 0) e B^, then Theorem l4.3l suflices for that value of x. Thus 
we only need consider values of x such that u{x, 0) > 0. For these values of x, we have 



which we choose to rewrite as 



\imy"dyU(x,y) = yu'^ ^ 



limy^dyiu^-^) = Cy 
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A bit of calculation shows us that m^"^ has V • (y''V(u^~'^)) > 0. Let w = u^~^ along y = Q 
and y = I and also along n {0 < y < 1 }, and satisfy V ■ {y"Vw) = inside. Then the 
maximum principle tells us that w > u^'^ since w^"^. Since u has a nontrivial free boundary, 
Theorem 14.31 tells us that u, and hence w, is bounded along y = I, say, w|y=i < C. Then 
that tells us that 

w{x,y)<w{x,Q) + C'y^'' 

inside ^3 n {0 < y < 1}. Hence, 

\u^-\x,Q) - u^-''{x,y)\ < Cy^" 
Applying an elementary inequality, we have the desired result. □ 

With these two lemmata in hand, we can prove the analogue of Theorem 14.121 for the 
domain where y > 0. 

Theorem 4.15. Suppose /3 < 1. Then there exists a K = K{6,n,j3), such that if Xi = 
(xi,yi),X2 = (X2,y2) 6 18-+^^ are in a 6-neighborhood of the free boundary, we have 

K^i) - m(^2)I < i^l^i -Xjf 

If IS > 1, there exists a K = K(6,n,fi, a), such that if Xi,X2 e K.+^^ are in a 6-neighborhood 
of the free boundary, we have 

\u{X,) - u(X2)\ < K\X, -X2r 

for any a < I. 

Proof. Without loss of generality, assume that yi < y2- 

First, assume that /3 < I. Suppose that 3^2 ^ l^i - ^2!'"^ • Then, using the C^-regularity 
of u restricted toy = and the previous lemma, we write that 

\u(xuyi) - u(x2,y2)\ < \u(xi,yi) - m(^i,0)| + \u{x2,y2) - u{x2,0)\ + \u{xuO) - U{X2,G)\ 

< C{/,+/^) + C\x,-X2f 

< 2C\X^ - X2\'^P + C\Xi - X2f 

< C\X,-X2f 

On the other hand, if yi > \Xi - X2\^~^, then we use our pointwise gradient estimates 
and the special properties of this case to write that 

\u(X,)-u(X2)\ < \Vu(X)m-X2\ 
u{X) 

< ^\X,-X2\ 

yi 
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where X is some point on the line joining Xi and X2. 

If yi < \Xi - Z2p-^ and y2 > \Xi - Zsl^-^, then we consider: 

\u(X,) - u(X2)\ < \u(xuyi) - u(x2, |Xi - X2n + \u(x2, \X, - X2\'-^) - u(x2,y2)\ 

The first term is controlled by the first method above, and the second term is controlled by 
the second method. 

For the case when > 1, simply let replace |Xi - Z2p~^ in the preceding argument by 
|Xi-X2|i-'^. □ 

5 Non-degeneracy 

Our goal in this section is to prove that energy minimizers of 

J(u) = ^ J" y"\Vu\^dxdy + J u^dx 

possess the property they are non-degenerate, which is to say that near the free boundary, 
they grow away from 0, and do not stay small. To be precise, our final theorem is 

Theorem 5.1. Let be a point of the free boundary ofu, a minimizer of J{u). Then there 
exists a constant C > such 

sup u > Cr^ 

Br 

Our strategy for proving this theorem is first to show that at a fixed distance away from 
the free boundary, there is a point which attains the desired growth. 

Theorem 5.2. Let xq € T be a point such that d(xo, F(u)) = r, where F(u) is the free 
boundary. Then, there exists a universal constant T(n, cr, 7) > such that 

u(xi) > Tr^ 

where |xo - xi| < ^. 

Proof As is typical, we shall rely on the scaling property of energy minimizers, specif- 
ically, that on jpu(xo + AX) is still an energy minimizer. Hence, we can assume 
d(xo, F{u)) = 1 and we only need to show that there exists an xi with 

u(xi) > C\ 

for some Xi with |jco - X]\<\. 

The standard Green's identity applied to some test function > with support com- 
pactly contained within Bi{xq) tells us that 

f u{\imy''dy(l)) - (l){\imy''dyu)dx = - f mV • (y"V(f>)dxdy 

Jr JBi(uo)n{y>0] 

1 
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We notice that, along F n Bi(xo), we have m > 0, and hence 



\imy"dyU = yu^ ^ 



Thus, we attain the condition that 

^ yu^'^cpdx 



^ uY\my"dy(pdx + J u(x)V ■ {y"V(())dxdy 



Now let us suppose to the contrary that there is no constant ci, that is to say, for any 
e > there is a minimizer such that lal < s inside Bi(xq). Since d(xQ,F(u)) = 1, we 

2 

have d{B i{xo), F(u)) < |, and we apply optimal regularity to bound the interior term: on 
5i(xo) r\ [y > 0}, we have 

u{x) < C 

Putting these conditions together, and we get the argument that 

J ye^~^(pdx < J slimy" dycpdx + J ■ (y"V (p)dxdy 

for arbitrarily small s. Since the left hand side becomes very large and the right hand side 
goes to zero, we have a contradiction: u cannot be made uniformly arbitrarily small inside 
B\{xq) n r, and thus there exists a constant t such that m > r at some point, which we call 

Xi. □ 

Now we begin the proof of Theorem l5.1[ which is essentially identical to that given in 
|l5l, and reproduced here for completeness: 

Proof. The proof is divided into two steps. 

Step 1. Let M be a local minimizer in Bm such that 

• is a free boundary point, 

• 5i(ei,0) nr c {u > 0} nr, 

• M(ei,0) = T > where t is the constant from Theorem 15 .21 known to be bounded 
both from above and from below away from 0. 

We claim the existence of A > and M > universal, the latter being large, such that 

sup M > (1 + A)t 

Suppose not. This implies the existence of a sequence of energy minimizers for our prob- 
lem, {uk)keii, satisfying the three listed conditions, such that 



lim sup u = T 
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From our regularity theorems, the family (uk)t is equicontinuous, and may be assumed 
to converge uniformly on every compact subset of R."^^ to a function Uoo which satisfies 
limy^oy"dyUco > 0. Moreover, Mco(-,0) has a maximum at ei, thus it is constant from the 
maximum principle. Hence Uoo = T, a contradiction because is a free boundary point. 
Step 2. Assume that is a free boundary point. As in L8J, we construct inductively a 
sequence of points (x,„),„ e R", such that 

• u(x„,+i,0) > (1 + A)u(x,n,0) 

• If r,„ = d{Xm, [u = 0}) and x^ is a free boundary point realizing the distance, we 
have x,n+i e BMr„(x,„) with u(x,„+i,0) > Tr^„. This is from the construction of Step 1 
applied to the rescaling -jM(i"m + V, rmy). 

In particular, we have 

\Xm+\ - Xm\ < 2{M + \)r,n 

We end the induction at the first point Xn, which leaves B^. This is possible, since the 
sequence m(x,„, 0) grows geometrically in m, but is controlled by optimal regularity consid- 
erations. Let mo be the index of the first point to leave B^. Then we write 

niQ 

u{Xma+\, 0) = Y^{u{x,n+i, 0) - u{Xm, 0)) > /I ^ u{Xm, 0) 
m=0 



> 2j {" = 0} n 5i/ by Theorem 

— C ^ ^ ■^m^ 



> C" 

The last step is justified because C", A are both universal, and mo is bounded universally by 
the geometric growth of the construction. Hence, for all r > 0, we have 



sup 



u > C"y^ 



which by rescaling Mr to r was precisely what we set out to prove. □ 

Corollary 5.3. In terms of n-dimensional Lehesgue measure, the positivity set {u > 0} has 
positive density, bounded away from 0, in a neighborhood of any free boundary point. That 
is to say, 

\B, n > 0}| ^ ^ 

— > 6{n, 0-, r) > 

for any ball Br centered about a free boundary point. 

Proof. This is a consequence of nondegeneracy, which says that a sufficiently positive point 
exists, and of the Holder continuity of u (Theorem |4.12| ). □ 
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A Proof of Lemma 4.6 



That the Riesz potential of a radial function is radial is obvious from symmetry considera- 
tions. 

The proof of this lemma uses facts from the theory of Riesz potentials. The key facts 
we will use are given by the following theorem of Adams [|T]| (pp 772): 

Theorem A.l. /// e LP(W'), \ < p < oo, then 

1. la J e BMO if and only if M^J e L~(R") 

2. la J e BMO implies lai+aif ^ C"^ where < ori < 1 
where 



iMaJ)(x) = sup r"^'" \f(z)\dz 

r>0 Jb,.{x) 

is the fractional Hardy-Littlewood maximal function. 

The plan is to set ipix) = (1 - 'i\x\f'^"'xBi{^), show that tp e L'(R"), and subsequently 
that Me_2o-iA 6 L"", whence we can apply the theorem to get the desired result. 
First, we prove that ip e V: 

J i(f{x)dx = C„ Jjil-Sr-f-^^f^'-^dr 

< Cn f\l-3rf-^'^dr<C 
Jo 

since /3 -2cr > -cr > -1. 

Next, we consider the fractional maximal function. It is clear that the points of concern 
lie directly atop the singularity, that is, r = |. In a balls Bp about such a point, we see that 

r ifiix)dx < Cp"*^-^^ 

which is precisely the scaling needed to see that Mp_2(T^ ^ C. 

Hence, we can apply the theorem of Adams and we conclude our lemma. 
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B Some facts concerning degenerate elliptic equations with 
A2 weights 



These facts are classical for the case of uniformly elliptic divergence form equations (see, 
Li3jl)> and close analogues are apparent by following the work of Fabes, Kenig, and 
Serapioni lfT2ll . although they are not explicitly stated there. We demonstrate the connec- 
tion, and give precise statements, in the case of the equation 

V • (AVu) = 

where w(x)A\^\'^ < ^^A^ < w(x)A\^\'^ and w is some A2 weight. 

B.l The De Giorgi oscillation lemma 

Lemma B.l. Suppose u is a positive supersolution in B2 with 

\{x eBuu> 1}| > s\Bi\ 

Then there exists a constant C depending only on s,n, and cr such that 

inf M > C 

Although we prove the former statement, the obvious corollary concerning subsolutions 
follows from applying the lemma to 1 - m, and is what we actually use: 

Corollary B.2. Suppose u is a subsolution in B2 with 

\{xeBuu<0}\ >s\Bi\ 
and u < 1. Then there exists a constant < /i < 1 depending only on e,n, and cr such that 

sup u < p 

Bi 

2 

In line with fl^, we let w{B) = J^y'^'dxdy represent the integral of our weight over a 
ball. The proof of this lemma depends on a Poincare inequality: 

Lemma B.3. For any s > there exists a C{s, cr) such that for u e H^{B\) with 

\{xeBuu = Q}\ >s\Bi\ 

we have 

y"u^dxdy < C \ y''\Vu\^dxdy 

JBi 
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Proof. It is a classical result (see, for example, Kinderlehrer and Stampacchia, 11. A. 1 5 [T4]) 
that for smooth functions u on that vanish on a set of measure at least eBr, we have 



J Br \^ ~ 21 



However, this is precisely the point of departure for Fabes, Kenig, and Serapioni (Theorem 
1.2) [|T2l . where they prove for functions u satisfying this condition, we have 



I y"u^dxdy < Cr \ y"\^u\^dxdy 

%J By ^ By 



which is precisely the result we were looking for. □ 

The rest of the proof follows the proof given in [fT3l , which is reasonably short, so we 
reproduce it here. 

Proof. Assume a > 5 > - we will see that the final result is insensitive to 5, and so we 
can let 5 ^ 0+ at the end. 

Let V = (logw) , then v is a subsolution to the equation, bounded by log5 ^ Then we 
have (Theorem 2.3.1 in Oil) 

supv<c( — ^- — ( y^v^dxdy 

Bi \w(5i) Jgj 

Applying the Poincare inequality, we see that 



Bi 

1 



supv<C /IVvpJjcJj 



MB I 

We set the test function = ^ for ^ e €^(82). Then we obtain 

< f y"Vu ■ V{—)dydx = - f ^{Vuf + 2-Vu ■ V^dxdy 
J u J u 



whence we obtain 



By fixing ^ = 1 on 5i and giving it bounded first derivative, we have 



r y''\V{\ogu)\^dxdy <Cw{B2) 

JBi 



Combining our statements, we find 



sup(logM) < C 

Bi 

1 
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W(52) 

w(5i) 



1 



which gives 

mfu>e y^-'^'^'l 
1 

is bounded since all Ap weights have a doubling property (see [fT9ll . V.1.5), and hence 

inf M > C 

□ 



B.2 The De Giorgi-Nash-Moser Harnack Inequality 

Theorem B.4 (DeGiorgi-Nash-Moser Interior Harnack Inequality). Let ubea non-negative 
solution in Bi to the equation. Then for r < I, we have 

sup u < c(l - ry inf u 
where c, p do not depend on r or the center of the ball. 

Proof. This fact is a straightforward extension of the standard interior Harnack inequality 
(proved in [[T2|'). which simply states that, so long as the equation is satisfied in B2, we have 

sup M < C inf M 

1 - 

where C > 1 is invariant under translation or dilation of the ball. In what follows, we 
assume r > ^, since the standard inequality proves the result for the case r < \, and that 
the balls are closed. 

Suppose ^ > r > |. Consider the collection of balls Bi{x), where x e dB\_. The union 
of these balls, along with 5^(0), is precisely 5| . For every x e dB^, we have 

sup u < Cu(x) < C sup u 

B 1 (.r) B I 

5 2 

Let X* e Biht such that 5i(jc*) is a ball containing infg, u. Notice that 

2 4 J 

C sup u < C^u(x*) < inf u 



1 



3 



Hence, we have 



sup u < C inf u 

B, B3 
J 



We use the same argument to extend from to inductively, and we get 



1^- 



sup 



u < C^*^"' inf u 



Bt_2-k 



and so on, until we reach a the first k such that 1-2 > r. At this point, we recognize that 
k » - log(l - r). Plugging in, we get the desired result. □ 
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